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Abstract 

We compute complete spectra of the staggered lattice Dirac operator 
for quenched SU(3) gauge configurations below and above the critical 
temperature. The confined and the deconfined phase are characterized by 
a different response of the Dirac eigenvalues to a change of the fermionic 
boundary conditions. We analyze the role of the eigenvalues in recently 
developed spectral sums representing the Polyakov loop. We show that 
the Polyakov loop gets its main contributions from the UV end of the 
spectrum. 
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Introduction 

In recent years many numerical studies of low lying eigenvalues and eigen- 
vectors of lattice Dirac operators were published. The reason for analyzing low 
lying eigenvalues is twofold: Firstly, based on the Banks-Casher formula [1], a 
close connection of topological objects to the chiral condensate has been pro- 
posed [2j. Secondly, for low lying eigenvalues many results from random matrix 
theory are available and were tested numerically in detail. 

While the low lying eigenvalues and their relation to chiral symmetry break- 
ing are well analyzed, very little is known about the UV part of the spectrum and 
its possible connection to confinement. Partly this situation is due the fact that 
an evaluation of complete spectra of a lattice Dirac operator is a considerable 
numerical challenge. 

A new way to analyze a possible relation between Dirac eigenvalues and 
confinement was proposed in [3]. It was shown that the Polyakov loop can 
be written as a linear combination of spectral sums over moments of Dirac 
eigenvalues computed with different (fermionic) boundary conditions. For the 
quenched case the Polyakov loop P is an order parameter for confinement, 
with (P) = in the confined phase, while in the deconfined phase (above 
T c ) the Polyakov loop develops a non-vanishing expectation valued In |3J it 
was speculated, that the difference between confined and deconfined phase can 
be characterized by a different response of the Dirac eigenvalues to changing 
boundary conditions. This difference leads to a non-vanishing (P) only in the 
deconfined phase. 

In this letter we analyze whether this scenario can be established in a numer- 
ical study of complete spectra of the staggered lattice Dirac operator. Further- 
more we address the question which part of the spectrum, IR or UV, contributes 
most to the Polyakov loop. 

Spectral sums for the Polyakov loop 

The spectral sums for the Polyakov loop, first presented in [3], were derived 
for the Wilson Dirac operator. Here we use the staggered Dirac operator at 
vanishing quark mass, 
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where n and m are integer valued 4-vectors labeling the lattice sites and r/^in) 
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1 In the theory with dynamical fermions one could study correlators of Polyakov loops 
to analyze the static potential. 



1 



(— i) n i+---+»v-i is the staggered sign function. The U^(n) denote the SU(3)- 
valued gauge links and we have set the lattice spacing to a = 1. 

We use an L 3 x N lattice and require that the number ./V of lattice points 
in time direction (the 4-direction) is even. As it stands, the Dirac operator 
has periodic boundary conditions in all directions. Below we will also need 
temporal boundary conditions with a phase z or its complex conjugate phase 
z* . These are implemented by multiplying the last temporal links with z and z*, 
i.e., 1/4(71, n4 = N) — » z U4 (ft, 114 = N) or U4 (ft, 714 = N) — ► z* C4(n, ri4 = N). 
Here we make the particular choice of ^-valued boundary conditions and set 
z = e i27r / 3 . VVe consider the Polyakov loop averaged over all of space, 

1 N 

P = 7^E Tr 4 II U4(n,n 4 )} , (2) 

n 114=1 

where Tr c denotes the trace over the color indices. 

Following the arguments in [3J, one finds that the Polyakov loop is given by 
a linear combination of spectral sums, 

N . ^ / ,.,\ N , , / f4 \\ N~\ 
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Each sum runs over all 3L 3 N eigenvalues and AW,a£*\ and A^* denote the 
eigenvalues computed with periodic, z-valued, and z*-valued boundary condi- 
tions respectively. The Polyakov loop P thus is represented as a linear combi- 
nation of spectral sums for the iV-th power of the eigenvalues computed with 
three different fermionic boundary conditions in time direction. The boundary 
conditions for the gauge fields are always kept periodic. 

It is interesting to note, that exact zero modes do not contribute to the 
spectral sum for the Polyakov loop. Thus, isolated topological objects which 
give rise to a zero mode, do not play a role for building up the expectation value 
of the Polyakov loopH. 

We stress that the result © is an exact formula for the Polyakov loop of an 
arbitrary gauge configuration. Below we will study numerically the expectation 
value (P) of the Polyakov loop in a quenched ensemble. The r.h.s. of Q 
then is rewritten in terms of expectation values of the iV-th moments of Dirac 
eigenvalues computed with the three different boundary conditions, 
■>N 



(1) 



2 We remark that the staggered Dirac operator does not have exact zero modes, but 
for sufficiently smooth configurations the would-be zero modes can be identified relatively 
clearly and are very close to the origin [4]. 
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This formula relates the vacuum expectation value of the Polyakov loop, which 
originally is a purely gluonic quantity, to spectral sums of the Dirac eigenvalues. 
Since the Polyakov loop is an order parameter for confinement (in the quenched 
case), with (P) = in the confined phase (T < T c ) and (P) in the 
deconfined phase (T > T c ), the formula (HD allows to study the relation of 
confinement and the spectrum of the Dirac operator. 

As one crosses the critical temperature into the deconfined phase, (P) ac- 
quires a non-vanishing expectation value. Equation fl4|) implies that the response 
of the eigenvalues to the boundary conditions has to change at T c , such that 
the spectral sums on the right-hand side do no longer cancel and (P) / 0. In 
the subsequent sections we study the response of the eigenvalues to changing 
boundary conditions and explore which parts of the Dirac spectrum give the 
main contributions to the spectral sums. 

Distribution of the Dirac eigenvalues and their response to 
changing boundary conditions 

For a numerical study of the formula (j4j) we need to compute complete 
spectra of the Dirac operator using three different boundary conditions. Since 
a numerical evaluation of all eigenvalues is a demanding task we are restricted 
to relatively small lattices. Here we use quenched gauge field configurations 
generated with the Luscher-Weisz action [5j on lattices of size 6 3 x 4. 

We work with two different values of the inverse coupling, j3 = 7.6 and 
/3 = 8.0. Setting the scale with the Sommer parameter one finds lattice spacings 
of a = 0.194(4) fm and a = 0.135(1) fm for the two couplings |6]. With a 
temporal extension of N = 4 this corresponds to temperatures of T = 254 MeV 
for (3 = 7.6, and T = 364 MeV for j3 = 8.0. Thus, we have two ensembles with 
temperatures below and above the QCD phase transition which for the quenched 
case is at T c ~ 300 MeV. Using LAPACK routines we compute complete spectra 
of the staggered Dirac operator with the three boundary conditions for 2000 
configurations from each of the two ensembles. The statistical errors we quote 
for the averaged observables are evaluated with single elimination Jackknife. 

We begin the presentation of our numerical results with a discussion of 
the eigenvalue distribution. Since the massless staggered Dirac operator is an 
anti-hermitian matrix, it has eigenvalues on the imaginary axis. We order these 
with respect to their absolute value and the sign of the imaginary part (the 
eigenvalues come in complex conjugate pairs). For analyzing the distribution of 
the eigenvalues we divide |A| into small bins of size A|A| and count the number 
of eigenvalues, An, in each of the bins. In Fig.[T]we plot An/A|A| as a function 
of |A| for both ensembles below and above T c (see |7] for a similar analysis in 
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Figure 1: Distribution of the eigenvalues A as a function of |A| (for periodic 
b.c). In the l.h.s. plot we use MeV as unit for the eigenvalues, while on 
the r.h.s. the dimensionless quantity o|A| (lattice units) is plotted on the 
horizontal axis. 



the case of SU(2)). On the horizontal axis we use either physical units (l.h.s. 
plot), or the dimensionless combination a\X\ (r.h.s.). 

As can be seen from the r.h.s. plot, the curves for the distribution of the 
eigenvalues are very similar for the two ensembles when plotted in lattice units. 
At small |A| the density of eigenvalues is a little bit depleted for the ensemble 
with T > T c . This can be understood qualitatively from the Banks-Casher 
formula [I], which predicts a vanishing spectral density at the origin for the 
chirally symmetric phase (above T c ). The corresponding opening of a spectral 
gap is well documented in numerical studies |8J and this phenomenon is reflected 
in the lower density at small [A| seen in our data for T > T c . The area under the 
two curves has to be equal (= the total number of eigenvalues) and we observe 
a light enhancement of the density for the T > T c spectra near the maximum. 
When we use physical units on the horizontal scale (l.h.s. plot), the two densities 
are stretched with different factors due to the different lattice spacing for the 
two values of the coupling we use. 

It is obvious that the distribution of the eigenvalues plays an important 
role for the spectral representation ([4]). Energy ranges where the density of 
eigenvalues is large will in general be more important than those parts of the 
spectrum with a low density. In order to disentangle the role of the density from 
other aspects, such as the response to changing boundary conditions, we now 
study two observables for individual eigenvalues. 
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Figure 2: Average shift of the eigenvalues when changing the boundary 
condition of the Dirac operator, plotted as a function of |A|. In the l.h.s. 
plot physical units are used and lattice units on the r.h.s. 



The first property of the eigenvalues we consider is their average shift when 
the boundary conditions are changed. To quantify this effect we define the 

average shift s(A^), by comparing an eigenvalue for periodic boundary 

U) (i) 

conditions to its partners X z and A^*, computed with z- and z*-boundary 
conditions, 

S (A«) = (|A«-A«| + |A«-A»| + |A«-A«|)/3. (5) 

In Fig.[2]we show the average shift s(A) as a function of |A|, again using physical 
units on the l.h.s. plot and lattice units on the r.h.s. The most obvious feature 
of the plots is the fact that the shift of the eigenvalues is considerably stronger 
for the T > T c ensemble, in particular towards the IR end of the spectrum. This 
confirms an observation made in [9J, where it was shown, that for T > T c the 
size of the spectral gap strongly depends on the fermionic boundary condition. 
The plots furthermore show, that for T > T c the shift is stronger than for the 
data at T < T c , not only in the deep IR, but for all of the eigenvalues. 

Apart from the different total shift, the two ensembles display also common 
features: A clear maximum close to the IR end, a minimum for midrange values 
and another increase at the UV end. The plots demonstrate that the IR modes 
are shifted most, combined with a less pronounced shift of the largest eigenval- 
ues. It is interesting to note that the parts of the spectrum with large shifts 
coincide with low densities (compare Fig. [TJ). 

Finally, we remark that both curves show a drop for the first bin in the IR. 
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Figure 3: Contribution of an eigenvalue to the spectral sum for the Polyakov 
loop, as a function of |A| (l.h.s. plot: physical units, r.h.s. plot: lattice units). 



We attribute this to the would-be zero modes which are stabilized by topology 
and thus should not move at all. However, since the staggered Dirac opera- 
tor has only approximate zero modes, these eigenvalues move a little bit, but 
considerably less than the bulk modes, thus creating the drop in the lowest bin. 

Contributions to the Polyakov loop 

We have demonstrated that when changing the boundary conditions, dif- 
ferent parts of the spectrum are shifted by different amounts. However, in the 
formula (j4|) for the Polyakov loop the iV-th powers of the eigenvalues enter 
and the shifted spectra are weighted with the phases z and z*. Thus, we now 
consider the contribution c(X^) of an individual eigenvalue to the spectral sum, 



N 

(6) 



3iVL 3 

In Fig. Owe show the absolute value of the contribution c(A) normalized 
by the total Polyakov loopH as a function of |A|. For both ensembles the size 
of the contribution increases strongly towards the UV end of the spectrum. 
Furthermore, when plotted in lattice units (r.h.s. plot), the two curves show 
a similar behavior. They both start out with a modest slope which, after a 
small dip (or shoulder) near a|A| ~ 1.5, turns into a steeper ascent. Fig. [3] 



3 We remark, that on a finite lattice the Polyakov loop does not vanish exactly also 
below T c . However, this "microscopic value" vanishes in the thermodynamic limit. 
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Figure 4: Accumulated contribution to the spectral sum for the Polyakov 
loop, as a function of |A| (l.h.s. plot: physical units, r.h.s. plot: lattice units). 

demonstrates, that although according to the last plot the IR modes are shifted 
more, the UV modes have a larger contribution to the Polyakov loop. 

Having studied the contribution of the individual modes, we finally have to 
fold in also the distribution of the eigenvalues analyzed in Fig. [U To take the 
spectral density into account we compute the accumulation of the contributions 
c(A), by summing all contributions up to a given value of |A|, i.e., we analyze the 
cumulated quantity S|A'|<|A| c (^')- ' n Fig. S we plot the absolute value of this 
quantity, normalized with the total Polyakov loop, as a function of |A|. As was 
to be expected, the IR modes do not contribute a lot and most of the Polyakov 
loop is carried by the UV modes. What is somewhat surprising, is the fact 
that the cumulation is not a monotonically increasing function, but has a dip 
which is particularly pronounced for T > T c . This shows that the buildup of the 
Polyakov loop from the spectral sums is not a simple linear process. The r.h.s. 
plot demonstrates that, when plotted in lattice units, the curves for the two 
ensembles are rather similar, as could already be concluded from the similarity 
of the behavior of the corresponding curves in Figs. [1] and [3l We remark that at 
the UV end both curves reach 1.0 to machine precision, which is a good check 
that the exact formulas ([3]), §4$ were implemented correctly. 

Phase of the Polyakov loop above T c 

Above T c the phase of the Polyakov loop has values close to the angles of 
the group center, i.e., close to 0, 27r/3 and Air/3. This can be seen in the plot 
on the very right in Fig. [5] where we show the complete spectral sum for the 
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Figure 5: The Polyakov loop in the complex plane for 20 configurations 
above T c as reconstructed from the spectral sums with 28%, 56%, 78% and 
100% of the eigenvalues (left to right). Note that the two plots on the l.h.s. 
have a different scale. 



Polyakov loop for 20 gauge configurations. For the truncated sum, we observe 
a surprising phenomenon: When including less than roughly two thirds of the 
eigenvalues, one finds that the results for the Polyakov loop show a phase shift 
of 180 degrees (see the two plots on the l.h.s. of Fig. [5]). Only the largest 
third of the eigenvalues has the correct phase and, since the UV contributions 
dominate, determine the final phase (the two plots on the r.h.s. of Fig. [5]). 

In Fig. [6] we show the phase shift A(p of the truncated spectral sum relative 
to the phase of the full sum. It is obvious, that when truncating the sum at less 
than two thirds of the eigenvalues, the phase is shifted by a value of 180 degrees 
(i.e., a shift of tt), and only the UV eigenvalues drive the relative shift to zero. 




Figure 6: Phase shift of the truncated spectral sums for the Polyakov loop. 



8 



It is interesting to note, that the position where the phase starts to come out 
right coincides with the position of the dip observed in Fig. 01 To summarize, 
the truncated sum for every configuation first evolves in the opposite direction 
in the complex plane and after two thirds turns back to approach the correct 
full Polyakov loop. 

Summary and interpretation 

We have generalized the discussion of spectral sums of Dirac eigenvalues 
representing the Polyakov loop |3J to the case of the staggered Dirac operator. 
In order to study these spectral sums numerically we computed complete Dirac 
spectra with three different fermionic boundary conditions, using quenched en- 
sembles below and above the QCD phase transition. 

Different aspects of these spectra were studied, in particular the distribution 
of the eigenvalues and their shift under a change of boundary conditions were 
analyzed. Concerning this shift we established that the IR modes are shifted 
most and also towards the UV end found a small increase. A comparison with the 
results from the distribution analysis showed that the eigenvalues with largest 
shift coincide with the regions of lowest eigenvalue density. Qualitatively this 
pattern holds for both ensembles, but for T > T c the shift of the eigenvalues 
is considerably larger than for T < T c . This enhanced response to changing 
boundary conditions leads to a non-vanishing (P) for T < T c , while in the 
confined phase only a microscopic value of (P) emerges which vanishes in the 
thermodynamic limit. 

Concerning the buildup of the Polyakov loop, we considered the contribu- 
tion of an individual eigenvalue as well as the accumulated contribution. The 
contribution of individual eigenvalues does not take into account that the den- 
sity of the eigenvalues is a function of their size and thus disentangles the two 
effects of a varying density and the different contribution of individual eigen- 
values. Both the individual as well as the accumulated contributions show that 
mainly the eigenvalues in the UV build up the Polyakov loop. For the phase of 
the accumulated contribution we have shown that in the IR a phase shift of 180 
degrees is observed which vanishes at the UV end. 

We have established the following qualitative scenario for the relation be- 
tween Dirac eigenvalues and the Polyakov loop below and above T c : In both 
phases the eigenvalues respond to changing the boundary conditions, but the 
response is considerably larger in the deconfined phase. As a consequence, the 
linear combination of the spectral sums with different boundary conditions leads 
to (P) > in that case. Concerning the role of different eigenvalues, we find 
that the spectral sums for (P) are dominated by the UV end of the spectrum. 
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